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Abstract. Current disruption (CD) and the related kinetic 
instabilities in the near-Earth magneto sphere represent phys- 
ical mechanisms which can trigger multi-scale substorm ac- 
tivity including g lobal reorganizations of the magneto sphere. 
Lui et all (|2008|) proposed a CD scenario in which the ki- 
netic scale linear modes grow and reach the typical dipolar- 
ization scales through an inverse cascade. The experimen- 
tal verification of the inverse nonlinear cascade is based on 
wavelet analysis. In this paper the Hilbert-Huang transform 
is used which is suitable for nonlinear systems and allows 
to reconstruct the time-frequency representation of empiri- 
cal decom posed modes in an adaptive manner. It was found 
that, in the Lui et all (I2QQ8' 



► event, the modes evolve globally 
from high-frequencies to low-frequencies. However, there 
are also local frequency evolution trends oriented towards 
high-frequencies, indicating that the underlying processes in- 
volve multi-scale physics and non- stationary fluctuations for 
which the simple inverse cascade scenario is not correct. 



1 Introduction 

A key element of substorm physics is an enhanced trans- 
port of magnetic flux from the dayside magnetopause over 
the pole s to the magnetotail during the growth phase of sub- 
storms 



e. 2. Baumiohann et al. , [mi. In the tail the mag- 



netic energy is accumulated and subsequently abruptly re- 
leased during the expansion phase. In general, substorms 
are associated with multi-scale processes including electron- 
(tens of kms), ion- (hundreds of kms), MHD- (> thousands 
of kms) scales and global reconfigurations of t he ma gneto- 
sphere iNakamura et al. . . l2006l : \Laitinen et all l2Q07h . For 
this reason, the identification of location(s) and timing of 
substorm triggering mechanisms was not fully conclusive. 



Magnetic reconnection (MR), converting magnetic energy 
into kinetic energy of accelerated plasma and energetic parti- 
cles at 20 to 30 could explain many substorm signatures, 
such as magnetic field dipolarization due to flow breaking at 
the barrier of strong dipolar field (at ^ IORe), flux pile-up 
and the tailward motion of dipolarization front, the forma- 
tion of substorm current we dge, and the associated aurora l 



and ground based activities ( e. g. Baumiohann et ali Il999|) . 



On the other hand, the substorm triggering scenario called 
current disruption (CD) invokes kinetic instability processes 
leading to cross-tail CD (current reduction), dipolarization, 
formation of the substrom current wedge in the near -Earth 
region at ~ 10 Re and other auroral and ground-based sub- 
storm signatures Ql2004 . The CD scenario offers a pos- 
sibility to explain both the dipolarization and the tailward 
movement of dipolarization front without invoking recon- 
nection associated flow breakings or flux pile-up . In fact, 
some dipolari zations were obs erved without local Earthward 
plasma flows ( Lui et aliuOOw . Though the multi-spacecraft 
timing observations from THEMIS spacecraft ind icate pre- 
dominant substorm triggering due to MR (Ang elopouloA 
I2OO8), both mid-tail (MR) and near-Earth (CD) processes 
can contribute to the global substorm activity. For example, 
MR associated Ea rthward flows might trigger near-Earth CD 
ie.Q. Runov et ali hoO^). 

CD is associated with kinetic inst a bilitie s in rather local- 
ized regions. It was shown by ('2004') that a local ki- 
netic theory of instabilities can account for the excited waves 
exhibiting large growth for long enough time even within a 
thin current sheet. The stability analysis of cross-field ion- 
drift driven instabilities based on a two-fluid approach also 
showed that waves initially excited near the ion cyclotron 
frequency (~0.1 Hz, 10 s) can grow and reach the typi- 
cal time scales associated with dipo larization (^0.006-0.008 
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Hz, 120-170 s). It was suggested [Lui et ali l2008h that the 



corresponding experimental signature would be a develop- 
ing inverse cascade in wavelet time-frequency representa- 
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tion, i.e., the evolution of wave modes from higher (shorter) 
to lower (longer) frequencies (time scales) in time. Con- 
versely, forward cascades in turbulence transfer energy from 
large scales to small scales. The dynamics, driving sources 
and dissipation of turbulent fluctuations were reviewed by 
Borovskv and Funstem (l2003h . In the Earth's plasma sheet 
turbulence w as also observed wi thin reconnection associated 
bursty flows ( Voros et ~ai\ l2QQ4l) . However, developing tur- 
bulence was not identified, possibly because of the short du- 
ration of events or due to the observation of already fully 
developed turbulence with broad-band fluctuations (with no 
time evolution of characteristic frequencies). Observations 
of a developing turbulence would help us to identify the spa- 
tial or temporal scales of turbulent drivers (the scales of en- 
ergy input) and to recognize the direction of turbulent cas- 
cades (forward or inverse). Anyhow, this would require to 
observe the associated fluctuations over the characteristic 
scales of turbulent drivers first, followed by the development 
of multi- scale fluctuations over the inertial range of scales. 



2 Limitations of the wavelet analysis 

While the experimental identification of inverse cascades 
would be crucial for understanding the real growth rate of 
excited kinetic instability modes, the detection based on 
wavelet analysis has to be interpreted rather carefully. The 
main reason is that the widely used ba sic wavelet analy- 
sis is non-adaptive iliuang et al\, Il998h . Once the basic 
wavelet is chosen it is used for the whole data set. This 
is not a problem when observations with gradual frequency 
changes, chirp signals, or wave trains are analyzed. Ac- 
tually, the harmonic signals used to test the wavelet iden- 
tification of inverse cascade features by Lui et g/J ( 2008 ) 
belong to this category of well treatable signals. How- 
ever, sudden jumps (non-stationarity) or nonlinearity can 
generate spurious time-frequency-energy distributions in a 
wavelet re presentation. ^ The basic wavelet transform is 
also linear {Huang et (2/., 1998'). Nevertheless, wavelets can 
be used successfully for the multi-scale ana l ysis o f non- 
linear systems iKijewski-Correa and KareerrL bood) . The 
results of wavelet analysis, however, heavily depend on 
the mother wavelet, proper resolution scale discretiza- 
tion liKijewski-Correa and Kareenl |2006|) . resulting in non- 
adaptiveness for multi-component non- stationary data. The 
effect of non-stationarity is demonstrated in Figure la, b. 
The test data consisting of a single sine wave with period 
1 s and a sudden jump (Figure la) is analyzed using the 
Morlet wavelet. The time- wavelet period-energy distribution 
(wavelet scalogram) in Figure lb is showing that the time lo- 
calization of the sine wave (between 4 and 5 s) is good only 
at wavelet periods near ~ 0.125 s. The time localization 
is rather poor at the wavelet period of Is, which is the ac- 
tual sine wave period. Using a different, optimized wavelet 
might lead to a slightly better localization of the sine wave. 
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Fig. 1. a. Test signal: sine wave + sudden jump; b. Wavelet analysis 
of the test signal; c. Hilbert-Huang analysis of the test signal. 




Fig. 2. Top subplot: magnetic field Bz component; Bottom sub- 
plots: empirical modes; emis, for i = 1,2,... 6. 



Certainly, a mother wavelet cannot be equally good for an 
unknown combination of periodic or stochastic components 
in a time series. 

Sudden jumps lead to broad-band energy distributions 
(Figure lb). Since for higher frequencies the basic wavelet is 
more localized, the sudden jump at Time=9s is well localized 
in time over small wavelet period ranges. The time localiza- 
tion worsens towards longer periods where the basic wavelet 
is more stretched. The non-stationarity also induces spurious 
periodicities in the time-period plane, in our case between 
Time=7-8 s, where the test signal contains no fluctuations at 
all. Actually, it can be even worse than that. The wavelet 
identification of scales (periods) can be entirely misleading. 
For example, Huang et al (1998) argue that a change occur- 
ring over large scales locally in time will appear in the cor- 
responding time-period distribution over small scales, due to 
the basic wavelet resolution. That is, a localized large-scale 
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change (non-stationarity) manifests itself over small scales in 
the wavelet spectrum. 

3 The Hilbert-Huang method 

To be able to identify inverse cascade features we need a 
method which is adaptive and suitable for non- stationary and 
nonlinear analysis of multi-scale data. The Hilbert transform 
(HT) represents a convolution of time series X{t) with time 
1/t, therefor e accentuates the lo cal properties of X{t). HT 



is defined as iHuang et ali Il998|) 



dr 



(1) 



where Pis the Cauchy principal value of the integral. X{t) 
and Y{t) are orthogonal and form a complex conjugate pair 



Z{t)=X{t)^iY{t) = a{t)e'^^'^ 



(2) 



a{t) and Q{t) are the instantenous amplitude and phase, 
respectively. The instantenous angular frequency uj{t) is the 
time derivative of &{t). HT is known for a long time, but 
the practical value of instantenous frequency uj{t) was ques- 
tioned, because of non-unique local values. Local unique- 
ness can be ach ieved by putting som e constraints on the data 
(Cohen, 1995). 'Huang et all {199$) proposed a practical de- 
composition algorithm, called empirical mode decomposi- 
tion (HMD), which removes the ambiguity in determination 
of uj. The extracted empirical modes (ems) fulfill certain 
conditions without any loss of essential nonlinear features of 
the original dat a and also ensure a daptiveness to local signal 
characteristics (Huang et <2/.Vl998). The HMD procedure fits 
splines to local maxima and minima of X{t), finds the mean 
Xm from the two spline functions and computes the residu- 
als Xr{t) = X{t) — Xm{t). These steps are iterated until a 
threshold is reached. The first mode is given by eml{t) = 
Xm{t). The next mode is found in the same way as eml, but 
replacing the original data X{t) with X{t) — eml{t). Suc- 
cessive emi modes are found by repeating the whole pro- 
cess, until the last residue contains only a simple trend. Even 
complex signals can be decomposed to a finite number of 
emis. For example, HT based on t he EMD was success fully 
applied to intermittent turbulence liHuang et ali l2QQ8l) just 
recently. The extracted em modes represents narrow band 
stationary Gaussian processes. Moreover, the modes are lo- 
cally symmetric relative to the envelopes defined by the local 
maxima and minima. This means that the local means are 
not computed over a time window or a local time scale, but 
are obtained directly from the local values of the envelopes. 

The instantenous frequency is determined by taking HT of 
each emi{t). The time and frequency localization is demon- 
strated for the test signal in Figure Ic. The test signal was de- 
composed into three emis (not shown). Contrary to wavelet 
results (Figure lb) the localization is rather good, without 
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Fig. 3. Power spectral density for magnetic field (red thick line) and 
empirical modes; the dashed red line shows the scaling region over 
the inverse cascade frequency range. 



spurious energy distributions in the time-frequency distribu- 
tion. Now, we are goi ng to apply the H ilbert-Huang method 
to the event studied hy lLui etal. ( 2Q08|) . proclaimed to be an 
example of inverse cascade observed by the wavelet method. 



4 Identification of an inverse cascade 

CD and the associated dipolarization was observed by the 
THEMIS A spacecraft at Xqsm ^ SRe, between 07:45 
and 07:50 UT on 29 January 2008. (In what follows, we 
will use decimal hours). Since the event occurred near a sub- 
storm expansion onset, CD is considered to be a substorm 
triggering agent. The inverse cascade was found to be mos t 
remarkable in the magnetic component Bz iLui et al.\^200S). 

The top panel in Figure 2 shows the Bz component of the 
magnetic field. Using the Hilbert-Huang method described 
above six emis were found which represent the decomposi- 
tion of Bz into 'monocomponent' modes. The modes are 
shown in the subsequent subplots in Figure 2. In Figure 3 the 
power spectral density (PSD) is depicted, computed for Bz 
(thick red line, for a better visibility shifted) and eml — 6. 
The monocomponent modes are narrow-band fluctuations, 
having finite bandwidth in PSD. For example, em2 exhibits 
the maximum power around the frequency ~4.510~^ Hz. 
Nevertheless, Figure 3 represents only a rough demonstra- 
tion of monocomponent modes. The modes in PSD corre- 
spond to a global description based on the Fourier transform. 
The monocomponent modes identified through the Hilbert- 
Huang method differ due to the localized constraints intro- 
duced in their determination explained above in Section 3. In 
spite of the finite bandwith the Hilbert-Huang modes can ex- 
hibit non-stati onarity through frequency and amplitude mod- 
ulation {yHuang et a/.lll998^) . 
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Fig. 4. Top: magnetic field time series; Middle: Wavelet time-period-energy representation; Bottom: Hilbert-Huang time-instantenous 
frequency representation; Thick dashed lines in the middle indicate the wavelet time-scale evolution in time; Line arrows show the locations 
of short time-scale (high-frequency) activations occurring simultaneously with the largest jumps in Bz- 



Following the \Lui et all (l20Q8k findings the inverse cas- 
cade should operate roughly over the frequency range from 
10~^ to 5.10~^ Hz. The magnetic field PSD shows scaling 
(thick red dashed line) over a narrower range from 3.10"^ to 
5.10"^ Hz. Though we do not intend to identify the inverse 
cascade on the basis of PSD, this scaling region could be 
associated with the developing nonlinear interactions. The 
power of eml — 3 modes over the scaling region is by 2- 
4 orders of magnitude smaller than the power of em4 — 6 
modes. Moreover, the em6 mode contains only a smooth 
trend, therefore, we will investigate em4 — 5 further. 

Figure 4 compares the wavelet time-period-energy spec- 
trum (middle subplot) and the Hilbert-Huang time-frequency 
spectrum (bottom subplot). Bz is in the top subplot. There 
are three vertical and inclined dashed lines in the wavelet 
plot. The lines are guides for eyes to see the frequency 
trend of activity in time. At least in three cases, there is 
an increased wavelet power which appears first at shorter 
time scales near the ion gyroperiod (~ 10 s). Power over 



longer time scales appears later. Our wavelet representa- 
tion of Bz is slightly different than the one published by 



Lui et all (|2QQ8|) . Nevertheless, the wavelet representation 
sh ows the s a me qu alitative features as the wavelet analysis 
in ^Lui et ( 2008 ) paper: increasing wavelet power bridg- 
ing gradually the scales between the ion gyroperiod (~ 10s) 
and typical dipolarization time (> 100 s). It is easy to no- 
tice that the enhanced small-scale wavelet power (indicated 
by dashed vertical lines in Figure 4) and the largest amplitude 
jumps in Bz occur simultaneously (vertical arrow lines). The 
dubious small-scale activity can occur due to a local sharp 
change (like in Figure 1), which might have nothing to do 
with the inverse cascade. Having in mind the limitations 
of the wavelet method, the uncertainty connected with sharp 
jumps in the data cannot be removed. Therefore, a wavelet 
independent test is needed. The bottom subplot in Figure 
4, corresponding to the Hilbert-Huang time-frequency spec- 
trum of em4 and em5 empirical modes (these are within the 
PSD scaling region, see Figure 3), seems to substantiate the 
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idea of multi-scale physics. We note, the emi spectra cannot 
be interpreted roughly before 07.74 and after 07.84 UT be- 
cause of interval finite- size effects. The first high-frequency 
activation (at ~ 3.10"^ Hz after 07.75 UT) in em4 occurs 
simultaneously with the beginning of small amplitude fluc- 
tuations in Bz- There is no increased power visible in the 
wavelet spectrum at aro und 07.75 UT in Figure 4, neither in 
the wavelet spectrum o uLui et all (l2008h . which could be as- 
sociated with those Bz fluctuations. In em4 the frequency 
decreases to 5.10"^ Hz at 07.78 UT then, after a short re- 
turn to - 10"^ Hz, it faUs down to - 3.10"^ Hz at 07.81 
UT. Overall, the instantenous em4 time scales change from 
30 to 330 s as time proceeds. During the interval of interest, 
emb varies between 8.10"^ and 10~^ (10-100 s). There are 
a few correlations between the wavelet and Hilbert-Huang 
spectra in Figure 4. For example, the enhanced small-scale 
activity in wavelet spectra is associated with high-frequency 
fluctuations in emb between 07.76 and 07.8 UT. The first 
two local peaks in emb seem to correlate with the first two 
small-scale activations in the wavelet spectra between 07.76 
and 07.77 UT. On the other hand, the local peak in emb be- 
tween 07.78 and 07.79 UT has no counterpart in the wavelet 
spectra. Moreover, there are low-frequency valleys in be- 
tween high-frequency local peaks, which are related to the 
occurrence of multi- scale interactions more complicated than 
a simple inverse cascade. In contrary to the wavelet method 
the Hilbert-Huang results indicate that the inverse cascade 
picture cannot fully explain the observed multi-scale fluctu- 
ations. This does not necessarily mean that the linkage be- 
tween local kinetic instabilities and large-scale processes is 
negligible. 



5 Conclusions 

In this paper we have shown that the typical limitations of 
wavelet methodology make the identification of CD associ- 
ated inverse cascades doubtful and misleading. The obser- 
vations of CD associated fluctuations are short in time, the 
magnetic field changes are highly dynamic. As it was out- 
lined in the paper, the Hilbert-Huang method is more suit- 
able for nonstationary and nonlinear processes. Its adaptive- 
ness allows to compute straightforwardly the instantenous 
frequency and reproduce the time-frequency representation 
of the time series based on decomposed empirical modes. 
In qualitative agreement with the wavelet results, the time 
evolution of instantenous frequency shows a general trend 
of raising magnetic fluctuations from high-frequencies to 
low-frequencies. However, the details of multi- scale fluctu- 
ations and the time locations of high-frequency activations 
are rather different in the wavelet and the Hilbert-Huang 
spectra. In the latter, there are also intervals where the fre- 
quency locally increases with time. In this paper we put a 
particular emphasis on the differences between the wavelet 
and Hilbert-Huang representations, demonstrating that for 



the CD associated non- stationary magnetic fluctuations the 
wavelet approach can lead to misleading conclusions. We 
believe the more adaptive Hilbert-Huang approach can re- 
construct the time evolution of instantenous frequency bet- 
ter, therefore, the identification of turbulent cascades occur- 
ring over the appropriate frequency ranges could be more 
straightforward. Since we were not able to identify a sim- 
ple inverse cascade from the data, the contribution of other 
physical processes to the observed fluctuations cannot be 
ruled out. For example, large-scale variations in the cur- 
vature of the ambient magnetic field lines driven by bal- 
looning instability prior to substorm associated dipolariza- 
tion onsets can le ad to fluctua t ions o ver the frequency range 
10"^ - 10"^ Hz iSaito et ali l2008h . This is also the fre- 



quency range over which the inverse cascade proposed by 



Lui et a l. (2008) should occur. Global instabilities, including 



magnetosphere-ionosphere coupling ( KaA l2007|) can also 
drive low-frequency waves. As a matter of fact, the iden- 
tification of low-frequency fluctuations is rather difficult or 
impossible when the whole length of event observations is 
comparable to the period of waves. It is even more diffi- 
cult to detect the non-stationarity of low-frequency compo- 
nents over short intervals using wavelets. In this paper it 
is shown that the low-frequency components of fluctuations 
(em4 — 5 modes in Figures 3 and 4 ) are non- s tation ary dur- 
ing the current disruption event of \Lui et all (l2008h . As it 
was mentioned in Section 2 non-stationarity of large-scale 
low-frequency waves manifests itself over high-frequencies 
in the wavelet spectrum, therefore the wavelet observations 
cannot be fully trusted. 

There exist also physical reasons which makes the identi- 
fication of turbulent cascades difficult. Inverse or forward 
cascades in turbulence are related to local mode interac- 
tions in the Fourier space. It is known that in the presence 
of a uniform magnetic field magnetohydrodyna mic turbu- 
lence is neither isotropic nor local in Fourier space iAlexakil 
[20071 There exist experimental indications that such nonlo- 
cal interactions can occur near plas ma boundaries with devel- 
oped gradients iVoros et al. J2OO7V The near-Earth location 
where CD and magnetic field dipolarization occurs repre- 
sents a dynamically changing boundary between dipolar and 
more stretched field lines. Spatial and temporal changes can 
be mixed up in one- spacecraft in-situ observations and the 
multi- scale signatures indicating the occurrence of an inverse 
cascade in the frequency domain misleading. 
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